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Schistosomes, unlike malaria parasites, are in their diploid stage when targeted by the human immune 
system. Diploids can be either homozygous or heterozygous. The difference has profound significance for 
developing immunity and yet has not previously been addressed. We examined the implications of zygosity 
on immunity to a diploid pathogen. Schistosoma japonicum and showed that the diploid state, and its 
associated heterozygous advantage, significantly affects the outcome of attack by the immune system and 
the accumulation of antigenic diversity in the parasite population. We demonstrate here that diploidy 
provides a novel means of immune evasion for diploid pathogens. 

Although both Schistosomes and malaria parasites have sexual cycles, the Schistosomes are in the diploid 
stages of their cycle when targeted by the human immune response, while malaria parasites are in their 
haploid stages. Since antigenic variation appears to be a major foil in the fight against infectious disease, it 
is important to consider the basic differences between haploid and diploid organisms as they relate both to 
immune evasion and to the accumulation of antigenic diversity in the population. 

Schistosomes are diploid pathogens that cause the second most deadly human parasitic disease after malaria. 
Vaccination, either alone or in combination with drug treatment, represents our best long-term hope for con- 
trolling schistosomiasis^ ^ but despite decades of research and development we still do not have an effective 
vaccine. Recent attempts to develop a vaccine against Schistosoma mansoni, using tegument proteins (Sm-TSP-2 
and Sm29) to stimulate an immune response, have met with some success^ ^ We are, however, still perplexed by 
the fact that vaccination with Sj-TSP-2 and Sj29, the orthologue to Sm-TSP-2 and Sm29^"^, kills only a small 
percentage of the parasites. One likely reason for the meager level of success appears to involve variations of the 
immunogens^'^. To date, most information about the antigenic variation of pathogens derives from the study of 
organisms whose life cycle is predominantly haploid^^"^^. One important question is whether our knowledge of 
the control of haploid organisms is, in fact, directly applicable to the development of diploid organisms. Diploids 
can be either homozygous or heterozygous. Scientists often do not agree about the relative advantages of the 
diploid and haploid states, although most do agree that the so-called "heterozygous advantage" benefits diploid 
organisms under certain circumstances (e.g. in the sickle cell genotype) and thereby affects the survival of the 
organism. 

Our study focuses on the fundamental biological impact of diploidy on the battle between the immune system 
and diploid pathogens. 

Results 

Expression of the Sj29 gene. We first isolated and sequenced the entire Sj29 gene. A comparison of its gene 
sequence with cDNA sequences indicated that the exon-intron structure of the Sj29 is identical to that of Sm29 
(Supplementary Fig. lA). We then established that both copies of the Sj29 genes were expressed in individual 
worms by sequencing cDNA from each of six. Four worms were heterozygous, containing two different Sj29 genes 
and transcripts (Supplementary Fig. IB). Only a single type of Sj29 gene and cDNA sequence was found in either 
of the other two worms, and it was concluded that they were homozygous for the gene (Supplementary Fig. IC). 
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Sequence polymorphism. We obtained total 268 full-length 
sequences of the Sj29 gene from 134 individual worms derived 
from a single field-isolate of S. japonicum cercariae. From these, 37 
different allelic forms were identified (Supplementary Table 1). Two 
types of analysis were done to establish the relationship among the 
alleles: neighbor- joining analysis (Mega 4) indicated that these could 
be resolved into distinct groups (Fig. lA); cladistic analysis 
(MacClade) indicated that most of the worms derived from a small 
number of lineages or clades (Fig. IB). Here we refer to three groups 
as clades A, B and M (a mixture of alleles that could not be 
definitively included in clades A or B). We will refer to the latter as 
"clade" M, although it is likely to represent more than one lineage, 
and avoid reference to it in situations where direct lineage is a factor 
in interpreting results. 

Immunization and challenge. We chose the most common allele 
from each clade for subsequent immunological analysis: clade A, 
Sj29-2; clade B, Sj29-4 and "clade" M, Sj29-20. These alleles 
represent 86.0%, 83.6% and 49.2% of their individual clade, 
respectively. The recombinant proteins were produced in Pichia 
pastoris for immunization experiments. Mice (n=10 per group) 
were immunized with the antigens formulated with adjuvant ISA 
70 M or the adjuvant alone and PBS as controls and challenged 
with 40 cercariae. Survival adult worms from each group were 
recovered six weeks following challenge and genotyped for Sj29 
gene. The surviving worms in individual groups were categorized 



as homozygotes (AA, BB) and heterozygotes (AB, AM, BM and 
MM). We will henceforth refer to worms carrying both copies of 
the Sj29 gene from the same lineage as homozygous. 

Vaccination with Sj29-2 provided significant protection against 
the A A genotype. The results showed that only 4 worms that were 
homozygous for the A clade survived in the vaccinated group, repre- 
senting 2.3% of total survivors, while 34 AA homozygotes survived in 
the adjuvant control group, representing 13.7% of total survivors 
(Fig. 2A, Supplementary Table 2). This indicated 88.2% reduction 
in AA worm burden. Vaccination with Sj29-2 provided no signifi- 
cant protection against heterozygous A genotypes with 33 surviving 
worms (37.1% of total worms). The experiment was repeated with 
like results (Fig. 2B). 

We then tried the same experiment using a Clade B antigen, Sj29-4. 
The number of Sj29-4 targeted BB homozygous worms was again 
significantly reduced (Fig. 2C). Hence, immunization with single pro- 
teins from either of the two major clades resulted in significant pro- 
tection from infection by homozygotes of their respective lineages ( AA 
or BB). On the other hand, exposure to either of the proteins had no 
effect on the protection of heterozygote infection. We also exposed 
mice to the Sj29-20 protein "clade" M. The results were consistent with 
those shown for other clades, but the number of homozygotes in the 
field-collected challenge isolate was too small to lend statistical support 
for the idea of clade specific protection (Supplementary Fig. 2). 

We immunized mice with a mixture of Sj29-2, Sj29-4 and 
Sj29-20 antigens and challenged these mice with field- collected 
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Figure 1 | Analysis of Sj29 sequence. (A) Neighbor- joining tree of Sj29 alleles and Sm29 (AF029222) based on p, the proportion of amino acid difference 
was generated by using Mega 4. (B) Cladogram of Sj29. 
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Figure 2 | Comparasion of the distribution of homozygotes and heterozygotes in immunized and control groups. Mice ( 10 per group) were immunized 
with the antigens (Sj29-2, Sj29-4 or Sj29-mix) adjuvanted by ISA 70 M, the adjuvant alone and PBS, respectively and challenged with cercariae. The 
survival worms from individual group were pooled for genotyping of Sj29 gene. The survival worms from were categorized as homozygotes (Homo-A, 
Homo-B), heterozygotes (Heter-A, Heter-B and AB) and others(Non-A, Non-B). The number of the surviving worms with various categorized genotypes 
in the immunization groups were compared to that in the control groups exception of the Sj29-mix group (D) where the data was converted before 
statistic analysis (see "Statistic Anaysis" in Method section). (A) Sj29-2 immunized group of the first experiment. (B) Sj29-2 immunized group of the 
second experiment. (C) Sj29-4 immunized group of the second experiment. (D) Sj29 multiple antigens immunized group. 



parasites. This greatly enhanced the vaccine's protective effect; it 
resulted in a 41.4% reduction in the overall worm burden in 
comparison to 26.2%, 23.6% and 14.8% reduction observed when 
using individual antigen for immunization (Supplementary Fig. 
3A). We then determined the Sj29 genotype of surviving worms 
(Supplementary Table 2). Significantly, immunization with the 
mixture significantly affected the survival of AB heterozygotes, 
along with the homozygotes (AA) (Fig. 2D). The presence of 
antibody to both of the expressed alleles appears connected with 
protection (Supplementary Fig. 4). Moreover, immunization with 



a mixture of antigens reduced the number of eggs in the liver by 
48.8% (Supplementary Fig. 3B). 

Infection by a clonal isolate. To control the genotypical complexity 
of the infectious agent, we cloned parasite lines from the field- 
collected isolate, as shown in Fig. 3A. Snails were infected with 
single miracidia. Then cercariae from snails were paired to infect 
rabbits and adult worms from each rabbit were genotyped. One 
pair showed only 2 alleles (i.e. Sj29-4/Sj29-4 homozygous female 
worms and Sj29-2/Sj29-4 heterozygous male worms). This clonal 
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Figure 3 | Establishment and infection of clonal parasites. (A) Establishment of clonal parasites. Snails were infected with single miracidia. Then 
cercariae from snails were paired to infect rabbits and adult worms from each rabbit were genotyped. One pair showed only 2 alleles (i.e. Sj29-4/Sj29-4 
homozygous female worms and Sj29-2/Sj29-4 heterozygous male worms). This clonal parasite of S. japonicum was designated Sj6M8. Eggs from this line 
were collected to establish clonal parasites. (B) Comparasion of the distribution of homozygotes and heterozygotes in Sj29-4 immunized and control 
groups infected with clonal parasites. N: the number of sequenced survival worms. 
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parasite of S. japonicum was designated Sj6M8. Eggs from this line 
were collected to establish clonal parasites. Mice were immunized 
with Sj29-4 recombinant protein (clade B) and then challenged with 
the clonal line Sj6M8. This resulted in a 43.1% reduction in worm 
burden and a 53.6% reduction of eggs in the liver (Supplementary 
Fig. 5). Of the survivors in Sj29-4 immunized group, 82.7% were AB 
heterozygotes. In contrast, the ratio of BB homozygotes and AB 
heterozygotes in the control groups was approximately equal 
(Fig. 3B). Hence, exposure to the clade B antigen resulted only in 
significant protection against the BB homozygote. Importantly, the 
greater extent of overall protection can be attributed to the relatively 
greater representation of BB homozygotes in the clonal population. 
This, as well as data shown in Fig. 2, suggest that minor sequence 
differences seen within a single Sj29 clade are not the result of 
immune selection and do not subvert the protective response. 

Discussion 

The response of diploid organisms to selective pressure is fun- 
damentally different from that of haploids. This difference is evident 
in the response of the immune system to the tegument protein Sj29 
and the progression of parasitic infection in immunized animals. 
Here we identified 37 different Sj29 alleles among 268 sequences. 
The ratio of non-synonymous to synonymous substitutions (dN/dS) 
in these variant forms, however, shows no indication of selective 
pressure on the protein sequence (Supplementary Table 3). 
Consistent with this result, the alignment of the Sj29 sequences indi- 
cated that the vast majority of sequences fit into distinct lineages with 
a low degree of homoplasy (sequence similarities that result from 
convergence rather than ancestral relationships). This would further 
indicate a lack of immune selection. In contrast, investigations of the 
diversity of blood stage P. falciparum antigens indicated both select- 
ive pressure and a high degree of homoplasy^^'^^. Do these observa- 
tions indicate a flaw in the algorithms meant to identify selective 
pressure, or that diploid organisms can evade immune response in 
ways that these formulas do not detect? 

We selected one representative protein from each of three lineages 
and studied the immunological response they elicited. Immunization 
with a single protein provided significant protection only against the 
development of Schistosomes that were homozygous for the same 
allele. A greater number of Schistosomes were killed when immu- 
nized mice were challenged with a clonal-isolate than a field-isolate. 
The percentage of clade -specific homozygotes killed, however, was 
the same regardless of the source of challenge. This addresses two 
issues: 1) the fact that same percentage of homozygous worms evaded 
the immune response regardless of their numbers in the challenge 
population indicated that their survival did not stem from an inad- 
equate, adaptive-immune response. Solely focusing on a means to 
increase the potency of the adaptive immune response (e.g., using 
different adjuvants) may not, therefore, lead to a more effective res- 
ponse; other factors appear to be involved in parasite survival. 2) The 
heterozygous worms have an advantage over their homozygous 
counterparts in certain natural situations. This advantage plays a role 
in preserving genetic diversity in the population. With regard to the 
latter, a predictive model suggests that although the percentage of 
each Sj29 clade within the population would rise and fall, none would 
become extinct (Supplementary Fig. 6). 

Immunization with a mixture of proteins representing three 
clades, followed by challenge with a field-isolate, supported the con- 
cept that products of both alleles need to be recognized by the 
immune system for a protective response to occur. In this experiment 
the survival of both homozygotes and heterozygotes was affected. 
The protection elicited by the mixture of all three proteins 
approached the sum of the protections offered by them individually. 
This indicates that the separation of immunogenic alleles according 
to their clade has immunological relevance. Clearly, a random 



selection of antigens for immunization would not have had an addit- 
ive effect on protection. 

To date, the accepted approach to preventing parasitic infection by 
means of stimulating the immune response does not take into con- 
sideration whether the pathogen under assault is haploid or diploid. 
This study is the first to emphasize the importance of zygosity on the 
immune response to diploid parasites. 

Methods 

Parasites and animals. A field-collected isolate of S. japonicum from Guichi County, 
Anhui Province, China was used in all the experiments. Female 6-8 week old BALB/c 
mice were obtained from the Songjiang Animal Facility of the Chinese Academy of 
Sciences of Shanghai. To establish clonal parasites, snails were individually exposed to 
one miracidia. Single snail-derived clonal cercariae were collected. Rabbits were 
infected with different cercaria pairs. The adult worms were genotyped and their eggs 
were used for establishing labaratory lines. All procedures performed on animals 
within this study were conducted in accordance with and by approval of the Internal 
Review Board of Tongji University School of Medicine. 

Allelic expression analysis. To estimate whether both the maternal and paternal 
alleles of Sj29 were expressed, single worm was cut into two parts for genomic DNA 
and total RNA extraction. Genomic DNA, cDNA and RNA (negative control) were 
used as PGR templates. The sixth exon with 180 bp was amplified by PGR. Sequence 
chromatogram traces of genomic DNA-amplified products were used to determine 
homozygote or heterozygote. The chromatogram traces of SNP site showed double 
peaks in heterozygote but single peak in homozygote. We compared the 
chromatogram traces of SNP site of cDNA with genomic DNA to estimate whether 
the maternal and paternal alleles of Sj29 were expressed. 

Characterization of SNPs in Sj29 open reading frame. We searched for SNPs in Sj29 
coding region by sequencing RT-PGR products after cloning into the pBlu 2KSM 
vector. The cDNA of 150 different individual worms were synthesized. The complete 
576 bp Sj29 ORF was amplified by RT-PGR. The resulting RT-PGR products were 
sequenced directly or sequenced after being inserted into vectors. Sequence 
chromatogram traces of RT-PGR-amplified products were used to determine 
homozygote or heterozygote. For every homozygote, three positive bacteria clones 
were sequenced, while for every heterozygote, three or more clones were sequenced 
until both maternal and paternal alleles were sequenced. The raw sequence data were 
aligned using Mega 4 software^^ DnaSP software^*^ was used to calculate n: an estimate 
of the average number of nucleotide substitutions between any 2 sequences; S: the 
number of polymorphic sites in the sample; and A: the number of alleles. The 
Neighbor-joining tree of Sj29 alleles and Sm29 (GenBank Accession No.AF029222) 
based on p, the proportion of amino acid difference was generated by Mega 4 
software. MacGlade 4.08 was used for cladistic analysis (http://macclade.org). 

Test of Neutrality. These rates of substitution within species were estimated using the 
method of Nei and Gojobori with the Jukes and Gantor correction^^ as implemented 
in the MEGA 4 software. The rates of nonsynonymous and synonymous substitutions 
were compared with the Z-test of selection (MEGA 4). 

Recombinant proteins preparation. Three majority alleles of Sj29 were selected for 
protein expression, Sj29-2, grouped in predominant population A; Sj29-4, grouped in 
predominant population B; and Sj29-20, grouped in predominant population G. The 
amino acid sequences of extracellular domain (Phe27-Lysl70) were back- translated 
into DNA sequences using codon usage frequencies for the yeast Pichia pastoris. The 
recombiant proteins were expressed and purified as reported previously^**. 

Vaccination and infection. Mice (n= 10 per group) were injected subcutaneously 
with 20 )j,g of single recombinant protein (Sj29-2, Sj29-4 or Sj29-20) or antigen 
mixture (6.7 [ig of each protein) in 30 [i\ PBS formulated with 70 [i\ ISA 70 M 
(SEPPIG, France) three times at 3-week intervals. PBS group and 70 M adjuvant 
group were set as control groups. On week 9, mice were challenged with cercariae of 
Guichi field isolate or a labaratory line Sj6M8. Each mouse was infected with 40 
cercariae. On week 15, adult worms were recovered by perfusion from the mesenteric 
vein and separated. The level of protection of vaccinated groups was calculated as a 
percentage based on the reduction in worm burdens or liver eggs per gram of liver, 
compared with those from the control group. 

ELISA. Serum samples were collected three weeks after the third immunization. 
Enzyme-linked immunosorbent assay (ELISA) was performed to detect the antibody 
levels of IgG against Sj29-2, Sj29-4 or Sj29-20. ELISA was performed as described 
before^**. Cut-off values were determined from the mean plus three standard 
deviations from the pre-immunization sera. 

Genotype of individual worm. We selected the trimorphic nucleotide site 162 to 
determine worm genotype, C162C for genotype AA, G162G for genotype BB, A162A 
for genotype MM, C162G for genotype AB, C162A for genotype AM and G162A for 
genotype BM (Supplementary Fig. 7). 
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Statistic Analysis. Data of worm burdens and liver egg burdens were expressed as the 
mean ± standard error. Multiple comparisons were analyzed by a LSD one-way 
ANOVA test, with ap-value < 0.05 considered significant. The genotype frequency of 
survival worms in different groups was calculated, respectively. For statistic analysis, 
it should be better to use Kruskal Wallis test. However, since the number of survival 
worms from individual mice was very small and quite variable from 2 to 33, we were 
not able to further categorize the survival worms with various genotypes and rank the 
data. Therefore, we pooled the survival worms of each group (10 mice) for Sj29 
genotyping and used chi-square test or Fisher's exact test for analyzing the ratio of 
various genotypes. To increase the validity of our conclusion, we repeated the 
experiments including use of the antigen with different alleles. For Sj29-mix 
immunized group (Fig. 2D), considering that the mixture of the antigens induce 
immunity that killed both homozygotes and heterozygotes, which led to significant 
reduction of total worms and equal-ratio reduction of worms with various genotypes, 
we converted the data before the statistic analysis. Firstly we calculated the percentage 
of worms with various genotypes in PBS control group using the formula: e.g. AB 
genotype worms% = AB genotype worms/total worms (i.e. 49/252 = 19.4%). Then we 
used the percentage to calculate the number of the infected cercariae with various 
genotypes. Since we used 400 cercariae for infection of 10 mice each group (40 
cercariae per mouse), the number of the infected AB genotype cercariae is 19.4% X 
400 = 78. Thus, of 78 infected cercariae, the survival AB genotype worms were 49, 41, 
and 15 in the PBS, adjuvant and Sj29-mix groups, respectively. The ratios of survival 
worms in various groups were further compared by chi-square test. 

Prediction Model. We established a predcition model to study the change of the 
frequencies of vaccine targeted alleles (clade A and B alleles) and untargeted alleles 
("clade" M alleles) during period of 10 schistosome transmission seasons. We 
consider clades A, B and M, with frequencies ai, bi and mi, respectively. According to 
Hardy- Weinberg equilibrium, the frequencies of six different genotypes AA, BB, 
MM, AB, AM and BM are a^ bj, m^ 2aibi, 2aimi and 2bimi. We consider the 
survival rates of vaccine targeted genotypes AA, BB and AB as a, (3 and y, respectively. 
In the survival population, the frequencies of clades A, B and M will be changed and 
can be calculated as the following formulas: 

ai = (a? ^ a + a^. ibj. i y + a^. imj. i )/(a^_ ^ a + bf _ ^ + m^_ ^ + 2ai. ibi. ^ y + 2ai_ im^ i + 2bi. im^ i ) ; 
bi= (bi_ 1 iS+ai.ibi.iY+bi.imi_i)/(af_ ^ a + b? ^ + mf_ ^ +2ai.ibi.iY+2ai.imi.i +2bi.imi.i); 
mi= (mf_ ^ +ai.imi.i +bi.imi.i)/(af_ ^ a + b^. ^ + mf_ ^ +2ai.ibi.iY+2ai.imi.i +2bi.imi.i). 
In this study, ai=0.55, bi=0.35, mi=0.13, a=0.32, (3=0.58, y=0.37, and thus we 
predicted the variation tendency of a^ bi and mi. 
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